Scale interactions and scaling laws in rotating flows at moderate Rossby numbers and 

large Reynolds numbers 
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The effect of rotation is considered to become important when the Rossby number is sufficiently 
small, as is the case in many geophysical and astrophysical flows. Here we present direct numerical 
simulations to study the effect of rotation in flows with moderate Rossby numbers (down to Ro ~ 
0.1) but at Reynolds numbers large enough to observe the beginning of a turbulent scaling at 
scales smaller than the energy injection scale. We use coherent forcing at intermediate scales, 
leaving enough room in the spectral space for an inverse cascade of energy to also develop. We 
analyze the spectral behavior of the simulations, the shell-to-shell energy transfer, scaling laws, and 
intermittency, as well as the geometry of the structures in the flow. At late times, the direct transfer 
of energy at small scales is mediated by interactions with the largest scale in the system, the energy 
containing eddies with k± fa 1, where _L refers to wavevectors perpendicular the axis of rotation. 
The transfer between modes with wavevector parallel to the rotation is strongly quenched. The 
inverse cascade of energy at scales larger than the energy injection scale is non-local, and energy is 
transferred directly from small scales to the largest available scale. Also, as time evolves and the 
energy piles up at the large scales, the intermittency of the direct cascade of energy is preserved while 
corrections due to intermittency are found to be the same (within error bars) as in homogeneous 
turbulence. 



I. INTRODUCTION 

Strong rotation is present in many geophysical and as- 
trophysical flows. Its effect is considered to become im- 
portant when the Rossby number (the ratio of the con- 
vective to the Coriolis acceleration, or the ratio of the 
rotation period to the eddy turn over time) is sufficiently 
small. The large scales of atmospheric and oceanic flows 
for example arc affected by the rotation of the Earth. 
The Rossby number for mid-latitude synoptic scales in 
the atmosphere is i?o « 0.1 [sj- In the Sun, the typical 
Rossby number in the convective zone is Ro « 0.1 — 1 
[sot . Furthermore, the Reynolds number {Re, the ratio 
of the convective to the viscous acceleration) in these sys- 
tems is also very large, and the flows arc in a turbulent 
state. 

Besides rotation, stratiflcation is also important in 
the atmosphere, the ocean, and other geophysical flows. 
Many studies have considered solely the effect of rotation 
in a turbulent flow, as a first step to gain better under- 
standing of the fluid dynamics of geophysical systems. 
For rapid rotation (very small Rossby numbers), signifl- 
cant progress has been made by applying resonant wave 
theory 0, [13, [H, [13 and weak turbulence theory [2^ . 
In these approaches, the flow is considered as a superpo- 
sition of incrtial waves with a short period, and the evo- 
lution of the system for long times is derived considering 
the effect of resonant triad interactions. This explains 
successfully the observed enhanced transfer of energy 
from the small to the large scales [1^ , and sheds light on 
the mechanism that drives the flow to be quasi-two di- 
mensional at large scales [4ll . l43j . Energy in three dimen- 



sional modes is transferred by a subset of the resonant 
interactions to modes with smaller vertical wavcnumbcr. 
Spectral closures [3j[iB| give similar results. 

However, resonant wave theories are only valid when 
the rotation period is much shorter than the eddy 
turnover time at all scales. For large Reynolds numbers, 
small scales are excited with a characteristic timescale 
proportional to the eddy turnover time, that decreases 
as the scales become smaller. Therefore the approxima- 
tions made in such theories can break down at sufficiently 
small scales, provided that the Reynolds number is large 
enough for these scales to be excited. How the results 
of resonant wave theory extend to the case of only mod- 
erate Rossby numbers but very large Reynolds numbers 
is still unclear. Several phenomenological theories have 
been developed to consider the case with large Re (see 
e.g. [l^, [4^) leading to different results for the 

scaling of the energy spectrum. 

In numerical simulations, the study of rotating turbu- 
lent flows is constrained by the computational cost of 
properly resolving the inertial waves and the resonant 
triadic interactions, together with the cost of resolving 
the small scale fluctuations when Re is large. Inverse 
cascades were shown to develop and anisotropics to ap- 
pear in low resolution (32^ and 64'^ grid points) simu- 
lations [1, [3, [2^ , either solving the equations of motion 
directly or using a subgrid model. Small aspect ratio 
boxes were considered in [39l [4l| allowing for an increase 
in resolution. Simulations at higher resolution were done 
later by [3| studying in particular the behavior of the 
shell-to-shell energy transfer. Recently, simulations with 
large Reynolds number and small Rossby number were 
performed using 128^ grids and 8th-order hyperviscosity 
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[18l | , thus confirming the dominant role of resonant triads 
for rapid rotation at large Re, although the results also 
suggest that resonant wave theories can be valid only for 
a finite interval of time. All these simulations also give 
different results for the scaling of the energy spectrum at 
scales larger than the forcing scale; it was shown in [40| . 
using a truncated model, that this can be the result of 
how all the relevant timescales arc resolved. 

Here we study the effect of rotation in a turbulent flow 
in high resolution direct numerical simulations (up to 
512'^ grid points). Simulations at this resolution were 
also performed recently in [s^; in this case, energy was 
injected at the largest scale available and the focus was 
solely on the scaling of small-scale fluctuations, showing 
depletion of the energy cascade and reduced intermit- 
tency. Our main objective, on the other hand, is to study 
the statistical properties of the fluctuations in flows with 
moderate Rossby numbers (down to Ro w 0.1) but at 
Reynolds numbers large enough to observe the beginning 
of a turbulent scaling at scales smaller than the energy 
injection scale. To this end, we use coherent forcing at 
intermediate scales, leaving enough room in the spectral 
space for an inverse cascade of energy to develop when the 
Rossby number is small enough. We also use the largest 
value of the Reynolds number allowed by our grid to ob- 
serve a direct transfer of energy at small scales. After 
describing the simulations, we study its spectral behav- 
ior, the shcU-to-shcU energy transfer, scaling laws and 
intermittency, and finally the geometry of the structures 
in the flow. 



II. NUMERICAL SIMULATIONS 

We solve numerically the equations for an incompress- 
ible rotating fluid with constant mass density, 

^ + u} X u + 2n X u= -W + i^V'^u + F, (1) 
ot 

and 

V • u = 0, (2) 

where u is the velocity field, u; = V x u is the vortic- 
ity, P is the total pressure (modified by the centrifugal 
term) divided by the mass density, and v is the kine- 
matic viscosity. Here, F is an external force that drives 
the turbulence, and we chose the rotation axis to be in 
the z direction: $1 = Qz, with Q the rotation frequency. 

The mechanical forcing F is given by the Taylor- Green 
(TG) flow [H 

F — Fq [sin(A:oa;) cos(fcoy) cos{koz)x— 

— cos(A:oa;) sin(fcoj/) cos{kQz)y] , (3) 

where Fq is the forcing amplitude. Although the forc- 
ing injects energy directly only into the x and y compo- 
nents of the flow, in the absence of rotation {fl = 0) the 



TABLE I: Parameters used in the simulations. N is the linear 
grid resolution, ko the wavenumber used in the forcing, v the 
kinematic viscosity, the rotation rate, fmax the maximum 
number of turnover times computed; Re, Ro, and Ek are 
respectively the Reynolds, Rossby and Ekman numbers. 



Run 


N 


ko 


V 






imax 


Re 


Ro 


Ek 




Al 


256 


2 


2 X 10" 


-3 


0.08 


45 


900 


4.50 


5 X 10" 


-3 


A2 


256 


2 


2 X 10" 


-3 


0.40 


45 


900 


0.70 


8 X 10" 


-4 


A3 


256 


2 


2 X 10" 


-3 


0.80 


45 


900 


0.35 


4 X 10" 


-4 


A4 


256 


2 


2 X 10" 


-3 


1.60 


45 


900 


0.17 


2 X 10" 


-4 


A5 


256 


2 


2 X 10" 


-3 


3.20 


150 


900 


0.09 


1 X 10" 


-4 


A6 


256 


2 


2 X 10" 


-3 


8.00 


185 


900 


0.03 


3 X 10" 


-5 


Bl 


512 


4 


8 X 10" 


-4 


0.40 


17 


1100 


1.40 


1 X 10" 


-3 


B2 


512 


4 


8 X 10" 


-4 


1.60 


25 


1100 


0.35 


3 X 10" 


-4 


B3 


512 


4 


8 X 10" 


-4 


8.00 


40 


1100 


0.07 


6 X 10" 


-5 



resulting flow is fully three-dimensional because of pres- 
sure gradients that excite the z component of the velocity 
[35I . [43 1 . The resulting flow has no net helicity, although 
locally regions with strong positive and negative helicity 
develop. It is also worth noting that this forcing injects 
zero energy in the fc^ = modes, whose amplification 
observed in the strongly rotating cases is only due to a 
cascade process. 

Two sets of runs were done at resolutions of 256'^ (set 
A) and 512'^ grid points (set B). The parameters for all 
the runs are listed in Table U With Taylor-Green forc- 
ing, the spherical shell in Fourier space where energy is 
injected has wavenumber kp ~ A/3fcoj or equivalently, at 
a scale Lp = 2'K/kp. For the runs in set A, kp « 3.5, 
and for the runs in set B, kp w 6.9; as a result, there 
is more room in spectral space for an inverse cascade to 
take place in the B runs. 

All the runs in set A were started from a fluid at rest. 
At t = 0, the rotation and the external forcing were 
switched on, until reaching a turbulent steady state, or 
until an inverse cascade was well developed in the case 
of large rotation rates. The runs in set B were done as 
follows. Run Bl was started from a fluid at rest and after 
turning on the rotation and external forcing, the run was 
continued to reach a turbulent steady state. Runs B2 
and B3 were started from a snapshot of the velocity field 
from the steady state of run Bl, and both runs were con- 
tinued until a new steady state was reached, or an inverse 
cascade developed. This latter method proved useful in 
saving computing time, as no differences were observed 
when comparing the late time evolution of the runs in 
the two sets. In all simulations, a dissipative range was 
properly resolved, and the time step was much smaller 
than all the relevant timescales. 

We define the integral and Taylor scales of the flow 
respectively as 

fE(k)k-^dk 
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FIG. 1: Time history of the energy for set A; Al (sohd), A2 
(dash-dot), A3 (dash-triple dot), A4 (long dash), A5 (dash), 
and A6 (dot); the inset shows a detail of the evolution at early 
times when waves prevail. Note the large increase in energy 
as Ko decreases. 



and 



A = 27r 



/ E{k)dk 



1/2 



(5) 



where E{k) is the energy spectrum. Since for large CI an 
inverse cascade develops, these two scales are useful to 
describe the evolution of characteristic scales in the flow 
with time. However, to avoid a time dependence of the 
Reynolds and Rossby numbers, we define for each run 
the Reynolds number as 



Re 



and the Rossby number as 



LfU 



Ro 



U 
2Q.Lf 



We also define the Ekman number as 



Ek^^ 
Re 



2nLi' 



(6) 



(7) 



(8) 



The turnover time at the forcing scale is then defined as 
T = Lp/U where U = ('f*'^) is the r.m.s. velocity mea- 
sured in the turbulent steady state, or when the inverse 
cascade starts. The amplitude of the forcing Fq in the 
simulations is increased as is increased in order to have 
[/ « 1 in all the runs. 



III. TIME EVOLUTION AND SPECTRA 

Figure [1] shows the time history of the energy in the 
runs in set A. Runs A1-A4 show a similar evolution, but 
runs A5 and A6 evolve differently. As the Rossby num- 
ber decreases, a transient develops in which the total 
energy oscillates with a frequency that increases with Q. 



FIG. 2: Time history of the energy dissipation rate (labels as 
in Fig. [1} the inset shows the evolution at early times. 




FIG. 3: Spectral index q as a function of time in run A6, 
in the isotropic energy spectrum E{k) (solid), in the E{k±) 
spectrum (dot), and in the E{k\\) spectrum (dash). Note that 
the energy is dominated by the orthogonal modes. 



It takes a longer time for this transient to decay as 
increases, and then the energy increases suddenly and a 
turbulent regime develops. An inverse cascade of energy 
is observed in run A6 after t « 120. The increase in the 
energy observed after this time is also accompanied by a 
monotonous increase with time of the flow integral scale 
L. Even in the runs in set B, that are restarted from a 
pre-existing turbulent steady state, long runs are needed 
to reach another turbulent state after turning on the ro- 
tation. As an example, in run B3 it takes « 20 turnover 
times for the transient to decay and for an inverse cascade 
of energy to develop. 

The energy dissipation rate 21/ J uj'^/2dV as a function 
of time is shown in Fig. [21 As the Rossby number de- 
creases, the peak of the dissipation rate is reached at 
later times, and then it saturates. Note that during the 
early transient in runs A5 and A6, the dissipation is al- 
most negligible, while in the saturated state the mean 
dissipation rate decreases slowly with decreasing Rossby 
number. 

The shape of the energy spectrum evolves with time, 
specially after the transient as turbulence sets in, and 
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FIG. 4: Isotropic energy spectra at late times in runs Bl 
(solid, t ^ 16), B2 (dot, t « 24), and B3 (dash, t ^ 40) at low 
Ro. Two slopes are given as a reference. The inset shows the 
isotropic energy flux for the same runs. 



later again as the spectrum becomes dominated by the 
contribution from the largest scales when the Rossby 
number is small enough for an inverse cascade to develop. 
Figure [3] shows the time evolution of the spectral index a 
(the exponent in the region of the spectrum with k > kp 
that follows a power law ~ fc") in run A6. Three curves 
are shown, which correspond respectively to the spectral 
index computed in the isotropic energy spectrum E{k), 
in the perpendicular energy spectrum E{k±) (where k± 
denotes the wavevectors perpendicular to ft), and the 
parallel spectrum i?(fc||) (where fcy denotes the wavevec- 
tors parallel to fJ). 

Before t « 80, we cannot recognize a power law in the 
energy spectra. After t 80, the spectral indices in E{k) 
and E{k±) grow monotonically from a value of —7 until 
reaching a plateau with a w —3 at f w 110. The energy 
spectra E{k), E{k±), and E{k^^) show wide and steep 
power law behavior from i w 80 to i w 120. During this 
transient, the energy flux is almost zero, as can also be 
expected from the small value of the energy dissipation in 
run A6 before t « 110 (Fig. [J). The end of the transient 
a,tt~ 110 and the plateau in a correspond respectively to 
the increase in the energy and in the energy dissipation 
rate showed in Figs. [1] and [21 The spectral index in 
i?(/c||) also has a plateau with a w —4.5. However, as 
the inverse cascade sets in and the energy piles up at the 
largest available scale in the system, the spectral index 
changes again and seem to slowly evolve towards a ~ —2 
in both E{k) and E{k±_). 

Note that the inverse cascade only starts after « 10 
turnover times after the turbulent state is reached at i « 
110. This can be understood as follows. The energy 
spectrum observed before f sa 110 has almost no flux. 
Nonlinear transfer of energy is required for the flow to 
become two-dimensional under the effect of rotation [Tsl . 
SH, [13 J and the nonlinear transfer is negligible until t w 
110. Then, after a few turnover times, the flow undergoes 
a transition and the inverse cascade sets in. 

The early transient is only observed in the runs in set 
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FIG. 5: E{k±) and Il{k±) (inset) at late times in runs Bl, 
B2, and B3. Labels are as in Fig. [4] 
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FIG. 6: n(fcx)/(?t^) at late times in runs Bl, B2, and B3. 
Labels are as in Fig. [J] 



A; since the runs in set B are started from a turbulent 
steady state. However, after the transient the spectral 
evolution of the runs in set A and B is similar. Since 
runs in set B have more scale separation for an inverse 
cascade to develop when Ro is small enough, we focus 
now on this set of runs. We show in Fig. 2] the isotropic 
energy spectrum at late times in runs B1-B3. While runs 
Bl and B2 show no growth of energy at scales larger than 
the mechanical forcing, except for some backscattering 
with a ^ fc^ spectrum, run B3 at late times is dominated 
by the energy in the k ^ 1 shell. At scales smaller than 
the forcing scale, the spectrum of run B3 is steeper than 
that of runs Bl and B2, and compatible with a ~ k~^ 
scaling. The inset in Fig. 2] shows the isotropic energy 
flux in the same runs. Note that in run B3, the flux at 
scales larger than the forcing scale is negative and ap- 
proximately constant, indicating the development of an 
inverse cascade of energy for small Ro. At smaller scales, 
the energy flux is positive. We thus conclude that in ro- 
tating flows, both the direct and inverse energy cascades 
can cohabit. 

The energy spectrum E(k±) is shown in Fig. [51 to- 
gether with the energy flux Il{k±). The spectrum and 
flux arc similar to the isotropic ones (indicating most of 
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FIG. 7: £'(fc||) at late times in runs Bl, B2, and B3. Labels 
are as in Fig. |4l 



FIG. 8: Shell-to-shell transfer function T{Q, K) &t Q = 40 for 
runs Bl (solid), B2 (dot), and B3 (dash) at late times. 



the energy is in these modes), and n(A:^) confirms the de- 
velopment of an inverse cascade of energy in k± at scales 
larger than the forcing scale in run B3, and a direct cas- 
cade at smaller scales with a fcj^ scaling. Figure [6] 
shows the energy flux T\{k±^ normalized by the r.m.s. 
velocity in each run. Note that the increase of the flux 
observed in the inset of Fig. [5] is only due to the increase 
in the energy of the system as the inverse cascade piles 
up energy at the largest available scale. As Fig. [H] in- 
dicates, the actual transfer of energy is slowed down by 
the rotation, and run B3 shows a smaller normalized flux 
than the other two runs at scales smaller than the forcing 
scale. 

On the other hand, there is no clear scaling in the small 
scales in E{k\\), nor an inverse cascade at large scales (see 
Figure[7l). The E{k\\) spectrum in run B3 is steeper than 
the E{k±) spectrum, consistent with the results shown 
in Fig. [3] for run A6 at late times. Slopes ~ k~^^^ and 
~ k^ are shown in Fig. [7] only reference. 



IV. ENERGY TRANSFER 

In this section we study the scale interactions and en- 
ergy transfer in rotating turbulent flows. A study of the 
energy transfer in this context, albeit at lower resolution, 
was done before by We will focus on runs Bl, B2, 
and B3, that have enough scale separation for direct and 
inverse cascades to develop when Ro is small enough. 
Similar results were obtained in the analysis of the runs 
in set A. 

To investigate the transfer of energy among different 
scales we consider the shell filter decomposition of the 
velocity field. 



u(x) = ^Ui^(x), (9) 

K 

where K denotes a foliation of Fourier space in shells. 



that for our purposes can be taken as spheres 



Uif (x) = ^ Uke" 

if<|k|<i<'+l 



cylinders [l| 



K<\'k±\<K+l 



or planes [l| 



Ui^ii (x) 



A'<|k|| \<K+1 



ik-x 



(11) 



(12) 



Then, we can define the shell-to-shell transfer between 
these shells as 



T(Q, K) = - Juk{u- V)uq dx^ 



(13) 



This function expresses the transfer rate of energy lying 
in the shell Q to energy lying in the shell K. It sat- 
isfies the symmetry property T{Q,K) = —T{K,Q) 0], 
and the numbers labeling the shells Q and K can cor- 
respond to any of the foliations of Fourier space listed 
above [l|. In particular, we will study the cases T{Q, K), 
T{Q±, K±), and T{Q\\,K\\). The energy fluxes discussed 
in the previous section can be reobtained in terms of the 
shell-to-shell transfer function as 



n(fc) 



(14) 



K=0 Q 



where again the wavenumbers k, K, and Q can corre- 
spond to different foliations of Fourier space depending 
on the subindex. 

Note that for the definition of the shells a linear bin- 
ning is used. Alternatively, the shells can be defined 
by a logarithmic binning of spectral space with intervals 
(7"A'o, 7"''"^A'o] for some positive 7 > 1 and for integer n. 
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FIG. 9: Shell-to-shell transfer function T{Qj_,K±) at Qj 
40 for runs Bl, B2, and B3. Labels are as in Fig. |8] 
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FIG. 10: Shell-to-sheU transfer function r(Q||,_R'||) at = 
40 for runs Bl, B2, and B3. Labels are as in Fig. |S] Notice 
these transfers are roughly 5 times weaker than in the 1. case. 



However, logarithmic binning cannot distinguish transfer 
between linearly spaced neighbor shells (from the shell K 
to the shell K+1) from the transfer between logarithmic 
neighbor shells (from K to '^K). If the cascade is the 
result of interactions with the large-scale flow (e.g., with 
modes with wavenumber kp associated to the external 
forcing) , the energy in a shell K will be transferred to the 
shell K + kp. Logarithmic binning does not distinguish 
this transfer from the transfer due to local triadic inter- 
actions that transfer the energy from K to 7A'. For this 
reason we use linear binning, but we note that care needs 
to be taken when using the word "scale" that implies in 
general a logarithmic division of the spectral space. The 
transfer among logarithmic shells can be reconstructed 
at any time later by summing over the linearly spaced 
shells. 

Figure [5] shows the shell-to-shell transfer T{Q, K) at 
K = AO for runs Bl, B2, and B3 at late times. The neg- 
ative peak to the left indicates energy is transfered from 
these i^T-shells to the shell Q = 40, while the positive 
peak to the right indicates energy goes from the Q = 40 
shell to those A'-shclls. In runs Bl and B2 the shcU-to- 
shell transfer peaks at \Q — K\ « kp- This was observed 



before in simulations of homogeneous turbulence 0, [3^ , 
and indicates that the energy transfer is local (the en- 
ergy goes from a shell Q to a nearby shell although 
the step in the energy cascade is independent of that 
scale and related to the forcing scale). However, in run 
B3 the transfer strongly peaks at jQ — A'| w 1. The same 
effect is observed in T{Q^_^K±) shown in Fig. O This 
indicates that at late times in run B3, the direct trans- 
fer of energy at small scales is mediated by interactions 
with the largest scale in the system, the energy contain- 
ing eddies with fc^ « 1 (see Fig. [5]). As a result, the 
timescale associated with the direct cascade of energy in 
fcx increases (and its flux reduces, see Fig. [6]), since the 
energy is transfered in smaller steps in Fourier space than 
in the case of the Bl and B2 runs. 

The shell-to-shell transfer T{Q\\,K^) at Q\\ = 40 for 
the same runs is shown in Fig. [TOl The dependence with 
the Rossby number of this transfer function is less dras- 
tic. In all runs, the transfer function T{Q\\, K\\) peaks 
at IQii — A'li I w kp. Since there is no inverse cascade of 
energy in , the energy containing scale in this direc- 
tion does not change as the Rossby number is decreased, 
and neither does the position of the peaks in T{Q\\,K\\). 
However, note the drop in the amplitude of the trans- 
fer in run B3 for all shells except the ones satisfying 
|Q|| — -ft^ii I = kp. As a result, for small Rossby number 
the transfer of energy between shells with Q\\ and K\\ is 
quenched except for the direct interactions with the ex- 
ternal forcing. Most of the interactions responsible for 
the transfer of energy to small scales between different 

shells are then interactions with the forcing. 

Figure [Til shows the transfer functions T{Q^^,K±) and 
T((3||,A'||) in nms Bl and B3 for all values of K and 
Q up to 40. In all cases, the white and black bands 
near K « kp and Q kp indicate a small amount of 
energy injected by the external forcing that is directly 
transfered to all wavenumbers up to w 30. For K and 
Q larger than kp, the figures confirm the results of the 
direct cascade of energy presented in Figs. [SlandfTUl For 
wavevectors perpendicular to H, as the Rossby number 
is decreased the peaks in T{Q±, K±) move closer to the 
diagonal K± = [Figs. [TlTa) and (b)], indicating the 
direct cascade in the perpendicular direction takes place 
in smaller k-steps given by the largest scale of the system. 
For all wavevectors, the energy in the parallel direction 
[sec T(Q||,A'||) in Figs. [TlTc) and (d)] is transfered to 
smaller scales, and the cascade step does not depend on 
the Rossby number. However, all transfer except the 
transfer with |Q|| — A'|| | = fcj? is strongly quenched in run 
B3. 

The development of a non-local inverse transfer can be 
observed in Fig. [TlTb) for A'x < kj and Qj_ < kj. The 
transfer is inverse, since below the diagonal Qj_ = K \_ 
regions with negative (dark gray and black) T{Q±, K±) 
can be observed. This means that energy is taken from 
e.g., K± = 20 and transfered to shells with Q± < kp. 
The transfer is also non-local, since this inverse trans- 
ference takes place between disparate scales. The non- 
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a) T(Kx,Qx) Run B1 (b) T(Kx,Q,) Run B3 




10 20 30 40 10 20 30 40 



(c) T(K,„Q„) Run B1 (d) T(K,„Q„) Run B3 




10 20 30 40 10 20 30 40 

Kn K„ 



FIG. 11: Shell-to-shell energy transfer functions T{Q±, K±) (a,b) and r((5|| , ify) (c,d) at late times in runs Bl (a,c) and B3 
(b,d). Notice the quenching of the transfer in case (d), except for the interactions with the forcing scale. 




FIG. 12: Second order longitudinal structure function S2{£±) 
(where £± denotes increments were taken in the direction per- 
pendicular to n) for runs Bl (soUd), B2 (dot), and B3 (dash). 



magnetohydrodynamics (MHD) [1, Near the diago- 
nal Q± = Kj_ the transfer is more complex. The inverse 
transfer superposes with a (smaller in net amplitude) di- 
rect local transfer (dark spots below and near the diag- 
onal, and light spots above and near it, for K± and Q± 
smaller than kp). This small direct transfer of energy 
at large scales is the result of a reflection of energy at 
K ~ 1, and was also observed in studies of the inverse 
cascade of magnetic helicity in MHD [J . The reflection of 
energy in Fourier space when it reaches the largest scale 
in the box suggests that the late time evolution can be 
dependent on the boundary conditions, a property that 
was already observed in simulations of two dimensional 
turbulence 0, [l^, [l^, [H, [2^ . In our case, the simulations 
do not contain a large-scale dissipation mechanism (such 
as a hypo-viscosity), and therefore energy piles up at the 
largest available scale until its growth is stopped by the 
(small-scale) dissipation. 



local transfer of energy in rotating turbulence shares sim- 
ilarities with the inverse cascade of magnetic helicity in 
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FIG. 13: Scaling exponents ("p for the steady state of runs Bl 
(+), B2 (*), and for run B3 at different times; t ^ 20 (0), t ^ 
25 ( X ) , f « 30 (□) , and t « 40 ( A) . The solid line corresponds 
to the scaling exponents given by the She-Leveque model [s^ , 
the dash line is the Kolmogorov prediction = p/3, and the 
dotted Hue is = p/2. 

V. SCALING LAWS AND INTERMITTENCY 

In this section, we consider the anisotropic inertial 
range scaling of the runs in Table [1] as described by the 
longitudinal velocity increments in the direction perpen- 
dicular to rotation, 

5u(x,£_l) =f • [u(x + ^f)-u(x)], (15) 

where f is a unit vector perpendicular to O. The lon- 
gitudinal structure functions Sp{£±) (with displacements 
along £±) can then be defined as 

Sp{i±) = {Su{^,£^)), (16) 

where the brackets denote spatial average. If the flow is 
self-similar, we expect S3{£±) ^ i±'^'', where are the 
sealing exponents. In isotropic and homogeneous hydro- 
dynamic turbulence, the Karman-Howarth theorem im- 
plies <S'3(^) ~ £, and the Kolmogorov energy spectrum fol- 
lows from the assumption 5*3 (^) ~ l^^^ [2l|. In practice, 
the spontaneous development of strong gradients in the 
small scales of a turbulent flow gives rise to corrections to 
this scaling, a phenomenon referred to as intermittency. 

From dimensional analysis, if the energy spectrum at 
small scales in rotating turbulence is E ~ kj^, we ex- 
pect S2 ^ i±- Figure [T2l shows the second order struc- 
ture function for runs Bl, B2, and B3 at late times out- 
side the wave regime when the turbulence has developed. 
At small scales for all runs, 52 ^ i'^, consistent with a 
smooth field in the dissipative range. At large scales, S2 
is larger for run B3 than for runs Bl and B2, a signature 
of the inverse cascade of eirergy and of the development 
of large scale structures in the flow. The scaling of runs 
Bl and B2 at intermediate scales is compatible with the 
Kolmogorov spectrum, while the scaling in run B3 is con- 
sistent with the ~ fcj^ energy spectrum. Note that such 
a scaling can be understood as a slow-down in the energy 




(5u /a 

FIG. 14; Pdf of the longitudinal velocity increments {£± = 377) 
for run B3 at different times; t ~ 20 (solid), t ~ 25 (dot), 
t ~ 30 (dash), and t ~ 40 (dash-dot); 77 is the Kolmogorov 
dissipative length. 



TABLE II; Characteristic scales and dimensionless numbers 
of the runs in set B. t is the time, L|| and L±_ axe the integral 
scales using respectively the -B(fc||) and E{k±) spectra, A is 
the isotropic Taylor scale, Ro\ is the micro-Rossby number. 



and n 


= 2C3 


-Cs. 










Run 


t 


ill 


Li_ 


A 


Rox 




Bl 


16 


1.5 


0.9 


0.29 


3.70 


0.23 ±0.01 


B2 


24 


0.9 


1.6 


0.31 


0.91 


0.24 ±0.01 


B3 


20 


2.6 


1.2 


0.50 


0.12 


0.19 ±0.02 


B3 


25 


2.4 


1.5 


0.55 


0.11 


0.26 ±0.02 


B3 


30 


2.1 


1.7 


0.59 


0.12 


0.26 ±0.05 


B3 


40 


1.9 


2.8 


0.53 


0.33 


0.24 ±0.02 



transfer rate because of interactions between waves and 
eddies (see e.g., [H,!!!, H^); such a slow-down is consis- 
tent with the results of the transfer function presented in 
the previous section. Considering the energy flux in the 
inertial range e is slowed down by (see e.g., [lH, H^) 

e-^ Suj^Tn/r^^, (17) 

where ~ and ti^ ~ £±/Sui^ is the turnover time 
of eddies in the plane perpendicular to f2; the scaling 

Suj^ ~ £± (18) 

follows. 

Figure [T51 shows the scaling exponents up to order 6 
computed in runs Bl, B2, and B3. The scaling exponents 
are defined as the exponents in 

Sp{e^)^e'^l (19) 

in the inertial range associated to the direct cascade of en- 
ergy (i.e., for £± < Lp). Runs Bl and B2 behave as non- 
rotating turbulence, with Kolmogorov scaling ((p « 2/3) 
and intermittency corrections (the prediction ~ p/3 
of Kolmogorov, and the intermittency model of intermit- 
tency in homogeneous and isotropic turbulence of She 
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FIG. 15: (Color online) Three dimensional rendering of the 
vorticity intensity in a subvolume of 256 x 512 x 512 grid 
points of run B2. The top view shows the subvolume in the 
direction of the axis of rotation; in the bottom view the red 
and blue arrows indicate respectively the x and z axis. 

and Leveque [11] are shown in Fig. [T3] as a reference). 
However, run B3 has a distinct behavior with C2 ~ 1- 
As time evolves in this run, and the energy piles up at 
^± ~ I7 the second order scaling exponent slowly con- 
verges to this value. Low order moments follow the curve 
Cp — p/2; but high order moments deviate from the 
straight line. The level of intermittency in the flow in 
all these runs can be measured in terms of = 2(^3 — Qq. 
This quantity, together with the integral scales of the flow 
(based on the parallel and perpendicular energy spectra) , 
the Taylor scale, and the micro-Rossby number (based on 
the Taylor scale of the flow) , 



are given in Table [TTl for the runs in set B at different 
times. 

It can be seen that at late times run B3 evolves 
towards an anisotropic state in the large scales, with 



L±/L\\ w 1.5. However, at small scales the flow seems 
more isotropic and at late times (t ~ 40) in this run 
Aj^/A|| w 0.8. The micro-Rossby number in runs Bl, B2, 
and B3 take different values in the range 0.11-3.7. How- 
ever, the value of /i is, within error bars, approximately 
the same for all the runs. As a result, the intermittency 
in the direct cascade of energy seems to be independent 
of the Rossby number Ro and the micro-Rossby number 
Rox. 

Finally, Figure[T4]shows the time evolution of the prob- 
ability density function (pdf ) of the longitudinal velocity 
increments in run B3. Increments in the direction per- 
pendicular to rj were computed, and the increment was 
taken equal to three times the Kolmogorov dissipation 
scale 77 in each run. The velocity increments in each run 
were normalized by their corresponding root mean square 
deviation a. In agreement with the level of intcrmitcncy 
observed in the scaling exponents, the pdfs show expo- 
nential tails indicating a larger than Gaussian probability 
of large gradients to occur in the small scales. The am- 
plitude of the tails of the pdfs as a function of 5u/a does 
not change significantly with time. Moreover, the root 
mean square deviation a of the velocity increments 5u 
increases with time. So if the pdfs are plotted versus 5u 
(instead of versus du/a), the pdfs actually become wider 
at later times. This effect can be understood considering 
that once the inverse cascade of energy sets in, the total 
energy in the flow as a function of time increases. 



VI. STRUCTURES 

The intermittency reported in the previous section in 
the scaling exponents and the pdfs of velocity increments 
indicates that even after the inverse cascade sets in, the 
flow develops strong velocity gradients in the small scales. 
In this section, we present visualizations of the flow and 
consider the structures that emerge. 

Figure [15] shows a three dimensional rendering of the 
vorticity intensity in half of the computational domain 
(256 X 512 X 512 grid points) at late times. The top view 
corresponds to the subvolume in the direction of the axis 
of rotation. Only regions with strong vorticity are shown. 
Note that the flow is anisotropic and quasi-2D, as it is 
clear from the top view. In the bottom view, the develop- 
ment in the flow of large scale column-like structures can 
be seen. However, the columns display small scale struc- 
tures with thin vortex filaments. These filaments seem 
to be ordered according to the large scale pattern. The 
presence of regions with strong vorticity even when the 
Rossby number is small enough for the inverse cascade 
of energy to develop can be expected from the results 
shown in Figs. [T51and[T?l linked to the intermittency of 
the flow. 

The local relative helicity u) ■ u/(|a;||u|) in the same 
subvolume is shown in Fig. 1161 Unlike in isotropic and 
homogeneous turbulence, regions of strong vorticity are 
not correlated with regions of strong relative helicity. The 
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FIG. 16: (Color online) Three dimensional rendering of the 
relative helicity intensity in a subvolume of 256 x 512 x 512 
grid points of run B2. Blue corresponds to positive helicity, 
and red to regions with negative helicity. Only regions with 
\u> -ul/dwllul) > 0.95 are shown, indicating strong quenching 
of nonlinearities. Viewing points are identical to Fig. [151 

net helicity over the entire box averages to zero, and lo- 
cally regions of positive and negative helicity seem to be 
more isotropic and homogeneous than the other quanti- 
ties studied. Note that the local relative helicity can be 
quite ubiquitously strong, indicative of a local quenching 
of nonlinear interactions. 



VII. CONCLUSIONS 

In this work, we presented results of the study of 
the turbulent scaling laws and energy transfer in direct 
numerical simulations of rotating flows in periodic do- 
mains. Spatial resolutions of 256"^ (set A) and of 512^ 
grid points (set B) were used, while moderate Rossby 
numbers (down to Ro » 0.1) and large Reynolds num- 
bers (up to Re w 1100) were considered, with enough 
scale separation to observe both a direct and an inverse 



cascade of energy when the rotation was strong enough. 
Runs in set A were started from a fluid at rest, while runs 
in set B were restarted from a previous state of homo- 
geneous turbulence. In the former case, for Ro sa 0.1, a 
long transient was found in which the energy dissipation 
is small, as well as the energy flux to smaller scales. Dur- 
ing this transient, the energy spectrum has a wide but 
steep spectrum, and its slope monotonously increases as 
a function of time. After turbulence sets in and the in- 
verse cascade of energy develops, the energy spectrum 
evolves towards a. E ^ kj^ scaling at scales smaller than 
the forcing scale. This late time evolution is observed in 
both sets of runs. 

At late times, the energy flux in runs A5, A6, and B3 
indicates an inverse cascade of energy in k± at scales 
larger than the forcing scale, and a direct cascade of 
energy at smaller scales. The net flux to small scales 
decreases as the Rossby number decreases, while the am- 
plitude of the flux to large scales increases. No inverse 
cascade is observed in fc|| . These cascades were confirmed 
by the study of the shell-to-shell energy transfer. The di- 
rect transfer of energy at scales smaller than the forcing 
is local, although in the runs with small Rossby number 
the transfer in k± is significantly slowed down. In this 
direction, the energy is transfercd between shells K± and 
Q± with small steps given by — K±\ w 1. As a re- 
sult, the direct transfer of energy in k± at small scales 
is mediated by interactions with the largest scale in the 
system, the energy containing eddies with k± 1. The 
timescale associated to the direct cascade in k± then in- 
creases, and its flux reduces. In fc|| the transfer is direct 
at all scales, and a larger component than in the case of 
non-rotating turbulence is due to interactions with the 
forcing scale. These results are in good agreement with 
phenomcnological derivations of the energy spectrum in 
rotating turbulence that consider a slow down in the en- 
ergy transfer rate because of interactions between waves 
and eddies [H, \^ . The non-local interactions also lead 
to the development of anisotropies in the flow [i^ . 

The inverse cascade of energy that develops at scales 
larger than the forcing scale in runs A5, A6, and B3 is 
non-local, in the sense that the transfer of energy as- 
sociated to this cascade takes place between disparate 
shells in Fourier space. At late times, the inverse trans- 
fer superposes with a (smaller in amplitude) direct local 
transfer of energy. This small direct transfer of energy at 
large scales is the result of a reflection a.t k± = 1, when 
the peak of energy reaches the largest scale in the box. 
Consequently, the late time evolution of simulations of 
rotating turbulence may depend on the boundary condi- 
tions used, a property already observed in simulations of 
two dimensional turbulence [a, [H, HI, HI] ■ 

The study of structure functions in the direct cascade 
range shows that the second order scaling exponent for 
increments perpendicular to the rotation in runs with 
small Ro is C2 ~ 1, in agreement with the energy spec- 
trum. Low order moments follow the curve = p/2 but 
high order moments deviate from this law, an indication 
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of intcrmittency. The level of intcrmittency in the direct 
cascade of energy, as measured by the exponent ^, is the 
same for runs with and without rotation. The sponta- 
neous formation of strong gradients in the small scales is 
further confirmed by pdfs of the velocity increments and 
by visualization of regions of strong vorticity in the flow. 

More separation of scales is needed to study the intcr- 
mittency in the inverse cascade of energy. Because of its 
relation to small scale gradients, intermittency is believed 
to be associated only with the forward cascade of energy. 
The intermittency phenomenon is not observed in the 
velocity field in two dimensional turbulence for which 
the conservation of vorticity leads to an inverse energy 
cascade to the large scales [HI ES], although intermit- 
tency in the vorticity (which cascades directly to small 
scales) is observed. It is unclear how the dual cascade of 
energy (towards small and large scales) in rotating tur- 
bulence affects the intcrmittency in the inverse cascade 
range. While intermittency is associated with small scale 
events, in many cases the strong events can affect the dy- 
namics of the large scales, specially in systems close to 



criticality; as an example, intcrmittency is a possible ex- 
planation for the occurrence of extended minima in solar 
activity [l3, [l^l ; it is also known to affect the transport 
of momentum in atmospheric surface layers (27j . 
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